Entanglement in extended Hubbard models and quantum phase transitions 
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The role of two-point and multipartite entanglement at quantum phase transitions (QPTs) in correlated elec- 
tron systems is investigated. We consider a bond-charge extended Hubbard model exactly solvable in one 
dimension which displays various QPTs, with two (qubit) as well as more (qudit) on-site degrees of freedom 
involved. The analysis is carried out by means of appropriate measures of bipartite/multipartite quantum corre- 
lations. It is found that all transitions ascribed to two-point correlations are characterized by an entanglement 
range which diverges at the transition points. The exponent coincides with that of the correlation length at the 
transitions. We introduce the correlation ratio, namely, the ratio of quantum mutual information and single- 
site entanglement. We show that at T = 0, it captures the relative role of two-point and multipartite quantum 
correlations at transition points, generalizing to qudit systems the entanglement ratio. Moreover, a finite value 
of quantum mutual information between infinitely distant sites is seen to quantify the presence of off-diagonal 
long-range order induced by multipartite entanglement. 
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I. INTRODUCTION 

Among the ground-state properties of many-body quantum 
systems, correlations are recognized as fundamental to the 
comprehension of the critical behavior displayed at quantum 
phase transitions (QPTs). 1 In this respect, in recent years a 
crucial role has been played by the notion of entanglement be- 
tween subsystems and the related measures developed within 
the field of quantum information theory. 2 

The large amount of results achieved^— — relies on the ob- 
servation that a QPT is, in general, characterized through non- 
analyticities of the density matrix of an appropriate subsys- 
tem. The latter is the starting point for the determination of 
any measure of entanglement, i.e., quantum correlation, ei- 
ther within the subsystem or between the subsystem and the 
remaining system. In the latter case, Von Neumann entropy 
S is always able to capture the presence of a QPT (of finite 
order) for an appropriate choice of the subsystem. Also, some 
general information about the type of transition (for instance, 
its order) can be gained by looking at the type of singularity in 
S. 5 Nevertheless, in order to construct a complete description 
of QPTs, one should provide a plethora of other features, such 
as, for instance, critical exponents. This requires the evalu- 
ation of more punctual measures of entanglement, like that 
of quantum correlations between two points or among more 
subsystems (multipartite correlations). Such measures have 
already been investigated in relation to QPTs for qubit sys- 
tems. In Ref. |H it was shown that concurrence (measuring 
two-point entanglement) scales with universal exponent for an 
XY model, whereas in Ref. 6 the n-tangle measure was used 
to detect the singular behavior of multipartite correlations at a 
QPT. 

A certain number of interesting results have been obtained 
as well for QPTs and entanglement in correlated electron 
systems^- — the latter being, in principle, characterized by a 
larger number of degrees of freedom per site (typically 4) with 
respect to qubit systems. This point makes it necessary to use 
measures of quantum correlations which are tailored also for 



qudit systems. Most of such measures are difficult to eval- 
uate whenever the subsystem is in a mixed state, since t hey 
require an often out of reach optimization process. In Ref. [lO|, 
a method was proposed to distinguish at a given QPT the con- 
tribution of two-point entanglement from that of multipartite 
quantum correlations without entering the above difficulties. 
The method is based on the comparative use of single-site Von 
Neumann entropy Si and quantum mutual information 2y; 
the latter being a measure of all (quantum and classical) cor- 
relations between two generic sites The method provides 
a simple recipe: whenever the two measures display the same 
type of singularity at a given transition, the latter ascribed to 
two-point correlations; on the contrary, if the singularity dis- 
played by Si is seen differently by 2y, the transition is to 
ascribe to multipartite correlations. 



In the present paper, we investigate the critical behavior of 
entanglement measures underlying the above classification for 
an extended Hubbard model at T = exactly solvable in one 
dimension. This is achieved by using appropriate measures 
of two-point and/or multipartite quantum correlations devel- 
oped recently in quantum information, with particular empha- 
sis on negativity, 21 concurrence, 20 and entanglement ration In 
particular, in order to generalize the latter to qudit systems, 
the correlation ratio is introduced. Also, the relation between 
critical exponents at the transitions and the scaling behavior of 
the entanglement measure at those critical points is studied. 



The paper is organized as follows. In Sec. [TTJ we introduce 
the model and its exact solution; we also derive the one- and 
two-site reduced density matrices. In Sec. [ill] we describe 
different measures of bipartite and/or multipartite correlations 
for qubit and qudit systems. In Sec. |IV] we present and dis- 
cuss the results obtained for the various measures at the differ- 
ent metal-insulator-superconducting transitions which charac- 
terize the model. Finally, in Sec. [V] we summarize our main 
conclusions. 
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II. THE BOND-CHARGE EXTENDED HUBBARD MODEL is the spin label, a denotes its opposite, n ja = c) J(y c ja is the 

spin-cr electron charge, and (i, j) stands for neighboring sites 
The model is described by the following Hamiltonian: on the chain; u and x (0 < x < 1) are the (dimensionless) 

on-site Coulomb repulsion and bond-charge interaction pa- 
H BC = - t 1 ~ x ( ni ° + n os)\ c \a c ^ — A* n ia rameters; fi is the chemical potential, and the corresponding 

<t >cr ia term allows for arbitrary filling. 




where c ia . and c ia are fermionic creation and annihilation 
operators on a one-dimensional chain of length L; a =T,I 

I 
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FIG. 1: Ground-state phase diagram of H. LEFT: n-u plane; empty 
circles stand for empty sites, slashed and full circles stand for singly 
and doubly occupied sites, respectively. RIGHT: fi-u plane. 



The model is considered here at x = 1, in which case the 
number of doubly occupied sites becomes a conserved quan- 
tity and the role of spin orientation becomes irrelevant to many 
aspects: for an open chain, any sequence of spins in the chain 
cannot be altered by the Hamiltonian, whereas for periodic 
boundary conditions, only the sequences of spins related by a 
cyclic permutation can be obtained. In particular, the ground 
state turns out to be degenerate with the fully polarized state, 
and the system behaves as if at each site i the local space had 
dimension Di = 3. In practice, both the Hamiltonian and 
the local vector space can be written in terms of the Hubbard- 
like projection operators X± = |a)j(/3|j, with local algebra 
X"^X^ S — 5f3 1 X aS and nonlocal (anti-)commutation rela- 
tions given by 



Xfxf = (-)(<*+PU-y+6)x] s xf (2) 



here a — 0, 1,2, |0)j = |vac)i is the local vacuum, 1 1 ) ^ = 
Xl°\0)i is the singly occupied state (with odd parity), and 
\2)i = Xf°\0)i is the doubly occupied state. More precisely, 
as far as the ground state is concerned, the model Hamiltonian 
in the one-dimensional case can be fruitfully written as 



H= ^(X^-XfX^+HxO+n^Xf 

i i 

- ^+^)J2(X1 1 + 2XD . (3) 



In this form, H provides the full spectrum of Hbc at ac = 1 
for open boundary conditions and its full ground-state phase 
diagram for both open and periodic boundary conditions. 
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A. Spectrum and ground-state phase diagram 

The physics of the system described by H is basically that 
of jV s = {^2,. X^ 1 ) spinless fermions which move in a back- 
ground of L — N s bosons, of which N d = Xf 2 ) are dou- 
bly occupied sites and the remaining are empty sites. Both 
N s and N d are conserved quantities, and determine the total 
number of electrons N = N s + 2N d . 

The situation may be understood in the formalism developed 
by Sutherland in Ref. 14. We can say that, apart from constant 
terms, H acts as a permutator of just two Sutherland species 
(SSs), the N s fermions, and the L — N s bosons. In practice, 
empty and doubly occupied states — though different as phys- 
ical species — belong to the same SS, since the off-diagonal 
part of the Hamiltonian cannot distinguish between them. It is 
only the constant term counting doubly occupied sites which 
depends on the actual value of N d . 

The eigenstates are easily worked outr^i^ and read 



\ip(N s ,N d ) >-- 



-o-io 

' A k N , 



(4) 



the result also holds at finite L if periodic boundary conditions 
are chosen in Eq. (01. Here 



(L-N s ~ N d )\/(L - N s )\N d \ 



1/2 



is a normalization factor; Xl° is the Fourier transform 
of the Hubbard projection operator Xj°, i.e., Xl° = 



V\ exp{ij^jk)Xj° . Moreover, — J2f =1 X 2a is also 

known as the eta operator, commuting with H; (rf) Ni creates 
N d pairs (or doubly occupied sites). 

The actual ground state \iPgs(N s , N d )) is chosen among 
the eigenstates in Eq. by requiring that N s and N d mini- 
mize the corresponding eigenvalue 



E(N a ,N d ) 



~L~) 



+ uN d -ii{N s + 2N d ) . 



In Fig. Q] we report on the left the ground-state phase di- 
agram in the n-u plane (with n = N / L average per-site fill- 
ing); on the right part, the same diagram in the fi-u plane. The 
phase diagram presents various QPTs driven by parameters u 
and \i (or n). TableUgives the range of parameters character- 
izing each phase in the thermodynamic limit (L, N a — > oo, 
with n a = N a /L finite). Each transition is characterized by 
a change in the number of on-site degrees of freedom (DOF) 
involved in the state. Phase IV has just one DOF per site since 
each site is singly occupied. It is an insulating phase, with 
charge gap A^ = /i + — /i_ = u — 4, where it + is 
the energy cost for adding (removing) one electron. Phases I 
and I' (which is the particle-hole counterpart of phase I) have 
two on-site DOF: singly occupied sites and empty or doubly 
occupied sites respectively. This holds for phase III as well, 
where only empty and doubly occupied sites appear. Phase II 
is the only phase in which all three on-site DOF are involved. 



Region of the phase diagram 



GS energy 



I: n s = n 

n d = 
V: n s — 2 — n 

rid = n — 1 



U > U c {n) = — 2 COS 7T71 — u/2 

it > u c (n) n = 2 cos nn + u/2 
^ = 



—2/n sin (nn) 
-2/-7T sin (nn) + u(n — 1) 



II: n s — 1/tt arccos (— u/4) u g [— 4, w c (n)] 

n d = l/2(7i - n s ) 

III: n s = u < -4 a ± = q=(2 + u/2) 
rid = n/2 

IV: n s = n = 1 u > 4 fj,± = =p(2 - u/2) 
n d = 



-2/n^l- (u/4) 2 + u/2[n - 1/tt arccos(-«/4)] 
un/2 




TABLE I: Ground-state sectors and corresponding energies. Here u c (n) = —4 cos(7rn). Note that the values limiting the range of u and/or jj, 
in each sector are the critical values for the transitions. 



Note that, as far as the relevant physics is concerned, this 
seems to be related to the number of on-site SSs character- 
izing the phase rather than the number of on-site DOF. In fact, 
phases I, P and II — which all have both the bosonic and the 
fermionic SSs — fall in the Tomonaga-Luttinger class, since 



neither spin nor charge gap are present, whereas phase III, 
though characterized by empty and doubly occupied states, 
has just the bosonic SS; it is again insulating, with charge gap 
Aj n = -u - 4. 
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Despite the above observation, phases I, F differ form phase 
II since only the latter is characterized by the occurrence of 
off-diagonal long-range order (ODLRO) and superconducting 
correlations (which also survive in phase III): 



lim (Xf°X™ r ) = n d (l - n d ) 



(5) 



Note that ODLRO — though not allowing real superconduct- 
ing order at x — 1 due to spin degeneracy, which implies the 
vanishing of spin gap— — is at the very root of superconduct- 
ing order, which occurs at x ^ lJ^- 

Before discussing the various transitions in terms of the 
behavior of entanglement measures, let us recall some fea- 
ture of each of them in terms of standard theory. First of all, 
since Nd and N s are both conserved quantities, the transitions 
should be originated from level crossing. Indeed, they also 
occur at finite L. Nevertheless, none of them is of first or- 
der, since it can be easily checked from Table U that the first 
derivative of Eqs is always smooth. Second, N s (1 — N s ) 
and N d (1 — N d ) can also be considered as order parameters 
for the transitions, since all of the QPTs occur in correspon- 
dence with the vanishing of one or both the above quantities. 
Moreover, all the three transitions I—* IV, II—* III and II— » 
IV correspond to the opening of an insulating phase, charac- 
terized by a charge gap linear in u (and in p). This implies 
that the product of the dynamical exponent z and the critical 
exponent v of the correlation length is 1 for all the three tran- 
sitions. Furthermore, since p = (d+ z)v — 1, with p exponent 
characterizing the first derivative of the free energy, we ob- 
tain that all of the three above transitions have v = 1/2 and 
2 = 2. The situation is less clear at the transition II— > I,F, 
since on the one hand no spin nor charge gap opens in both 
phases. On the other hand, it must be said that Nd, and con- 
sequently r\ pairs and ODLRO, vanish in correspondence with 
the transition. This is true as well for the pairing gap Ap, 
A P = E(N + 2) + E(N) - 2E(N + 1). Indeed it can be 
seen that Ap = in phase I, and Ap = u — u c (n) < in 
phase II, where u c (n) defines the critical line. 



where q a p = (tpGs\X"^\tpGs)^ while the two-site reduced 
density matrix reads 



Pij = Tl "L/{y} (PGS) 



E 

a, ,3,7,5=0, 1,2 



with q a fjjs = (ipGs\X°' f3 Xj 5 \ipGs}- Below we report the re- 
sults for pi and pij (the detailed derivation of the calculations 
for the dimer case can be found in AppendixlBl. 

When expressed in terms of the basis {|0},, |2),}, pi is 
diagonal in all the regions of the phase diagram: 



diag {1 - n s - n d ,n s ,n d } 



whereas 



with 



respect 



to 



the 



(6) 
basis 



{|00), |01>, 1 10>, |11), |12), |21), |02), |20), |22)}, Pij reads 
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(7) 



with c = n d /(l - n s ), P l3 = (1 - n s ) 2 - \Y l3 \ 2 , e = c/L, 
and 



ir,. 



1 sin(n s 7r|i — j\) 
L sin(f|i-j|) 



B. Reduced density matrices 

The evaluation of the measures of correlation described in 
the following sections requires the manipulation of the single- 
site and dimer reduced density matrices, pi and pij . These 
can be obtained from the system density matrix in the ground 
state; the latter being defined as usual by pes = IV'gs) (ipGS I- 
The reduced density matrix pi (pij) is then the trace of pes 
with respect to all the DOF except those of the site i (sites 
i, j). These matrices can be constructed in a simple way from 
the one- and two-point correlation functions, using the opera- 
tor expansion for the density matrix in terms of the Hubbard 
projectors (O; one, however, has to pay attention to the exist- 
ing graded structure of the fermionic algebra, 18 (see Appendix 
lAl . In particular, the one-site reduced density matrix can be 
written as 

Pi = Tr L/l (p GS ) = ^2 qapX? , 

q,/3=0,1,2 



III. MEASURES OF ENTANGLEMENT 

The theory of quantum information has provided the study 
of complex quantum phenomena, such as QPTs, of new and 
well-defined tools. In many cases, these tools have been used 
to describe the (critical) behavior of relevant many-body sys- 
tems. In general, the fact that a critical behavior can be spotted 
by appropriate measures of entanglement should not be a sur- 
prise from the point of view of the Landau theory, since any 
measure of entanglement can be expressed as a unique func- 
tional of first derivatives of the ground-state energy. 5 Never- 
theless, the use of more advanced tools could provide new 
interesting features difficult to extract from standard theory. 
For instance, in Ref. 10 we have described how, by using the 
appropriate measures of bipartite correlations, it is possible 
not only to fully describe the phase diagram of some model, 
but also to discriminate the role of two-point from multipartite 
entanglement at each of the QPTs the system undergoes. 
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A. Separating two-point from multipartite entanglement at 
QPTs 

In this section we briefly recall the method used in Ref. 
where we were interested in the existing correlations between 
(a) the single site i and the rest of the system, and (b) the 
generic site i and a generic site j ^ i. 

Since the full system is in a pure (ground) state, the amount 
of quantum correlations between a single site and the rest of 
the system is measured by the Von Neumann entropy of pf. 

D 

S t =S(pi) = -^2\j\og 2 X j , (8) 

where Xj , j = 1 , . . . , D, are the eigenvalues of the reduced 
density matrix pi . 

The total correlations (quantum and classical) between two 
sites i, j are captured by the quantum mutual information: 

l t] = S( Pl ) + S{ P] ) - S( Pt] ) , (9) 

where S(pij) is the two-site Von Neumann entropy, S(pij) — 
J2j=i \j 1°§2 Aj> an d Aj , j = 1, . . . , D , are the eigenvalues 

Of pij. 

To resume, we have the following situation. On the one 
hand, the single site is quantum correlated with the rest of the 
system in two possible way: via two-point correlations (Q2), 
when it is "individually" correlated with some/all the other 
sites, and via multipartite correlations (QS), when it is con- 
nected through n-point quantum correlations. On the other 
hand, the mutual information allows one to evaluate all the 
correlations connecting two sites; the latter can be of a quan- 
tum nature, the already mentioned Q2, and/or of a classical 
nature (C2). 

While the single-site entanglement Si is not able to distin- 
guish between Q2 and QS, the quantum mutual information 
Tij is not able to distinguish between Q2 and C2. Never- 
theless, in Ref. 10 we have shown that a comparison of the 
singular behavior of Si with that of ly allows one to discrim- 
inate whether a QPT is ascribed to Q2 or QS correlations. In 
fact, whenever the singular behavior exhibited by Si is due 
to Q2 correlations, the same type of singular behavior is nec- 
essarily displayed by as well, since it also contains Q2 
correlations. 



B. Measuring two-point entanglement 

The task of measuring quantum correlations between two 
given sites i and j has a simple solution when i and j are 
two-level systems (qubit) in terms of the concurrence. 20 Even 
when i and j are arbitrary qudit, the quantification of entan- 
glement can be carried on by means of the negativity^ the 
latter being a lower bound for concurrence^ 



1. Concurrence 

The concurrence was first introduced in Ref. and, for 
the case of two qubits, it is directly related to the entangle- 
ment of formation. In order to evaluate the concurrence one 
has to first manipulate the two-qubit density matrix p^ and 
find p~ij — PijCTy <jyp*j<T y 0) o y , where p*- is the element- 
wise complex conjugate of pij. The concurrence can then be 
written as 

Cij = C(pij) = max{0, Xx - A 2 - A 3 - A 4 } , (10) 

where the A^'s are the square roots of the eigenvalues of plj 
taken in decreasing order. 

2. Negativity 

Another measure of bipartite quantum correlations Q2 is 
the negativity, 

U ij =Af( Pij ) = (\\pJj\\ 1 -l)/(d-l), (11) 

where p^ is the partial transposition with respect to the sub- 
system (site) i applied on pij, and ||0||i = TrvO^O is 
the trace norm of the operator O. pjj can have negative 
eigenvalues p n and the negativity can also be expressed as 
N{pij) = | 2~2 n Pn\- Although the negativity is not a perfect 
measure of entanglement, 30 since it fails to signal the entan- 
glement in the subset of mixed states called partial positive 
transpose states, it gives important bounds for quantum infor- 
mation protocols, i.e., teleportation capacity and asymptotic 
distillability. 

One would reasonably expect the measures of two-point en- 
tanglement to exhibit the same singular behavior of Si and 
when the transitions are ascribed to Q2 correlations. This is 
not always the case. 10 To understand such an unexpected fea- 
ture, we shall explore in more detail the behavior of both A/y 
and Cij when r = \i — j\ is varied in proximity of the QPTs 
dominated by Q2 correlations. 

C. Multipartite entanglement measurements 

The case in which the singular behavior of Si is ascribed to 
QS correlations can also be treated with the described bipar- 
tite measures in a simple, though not complete way. The only 
thing that one can say is that when the QS correlations enter 
into play, the same singular behavior should not be displayed 
by or by A/y, since both measures regard only two-point 
correlations. 

We now proceed to review the measures of multipartite en- 
tanglement useful for our analysis. 

1. Residual entanglement: The tangle 

The idea of residual entanglement was first introduced in 
Ref. where the case of a three-qubit system in a pure state 
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\^abc) was studied. The basic ideas are as follows: 

(i) The concurrence for a two-qubit pure state reduces to 

C A . B = 2 V 'det p A . (12) 

(ii) Once a focus qubit is chosen, in this case A, the follow- 
ing inequality holds for a three-qubit pure state: 

4det p A > C\ B + C\ c . (13) 

(iii) In the case of a three qubit pure state, the subsystem 
constituted by the pair (B, C) is four dimensional, but only 
two of these dimensions can be used to express the state; in 
other words, both the reduced density matrices pa and psc 
have only two nonzero eigenvalues. This fact leads one to 
interpret 2v / defpX as the concurrence between A and (B, C) 
and thus to rewrite the above inequality as 

C%{bc) ^ Cab + C\c ■ (14) 

The last result says that the entanglement that the focus qubit 
A can establish with each of the other qubits separately is 
bounded by the entanglement that it can globally establish 
with them. (The latter is a property that is not satisfied by 
the entanglement of formation.) The definition of residual en- 
tanglement for a three-qubit pure state, or tangle, can be intro- 
duced in the following unique way on the basis of the above 
results and of the fact that they do not depend on the focus 
qubit chosen: 

TABC = C 2 a ( BC) - C\ B - C\ c . (15) 

Due to the permutation invariance, this quantity properly mea- 
sures at least an aspect of three-qubit entanglement: the three- 
way entanglement. 

2. Entanglement ratio 

An example of the use of the tangles for the exploration 
of QPTs in spin systems is given in Ref. 0. There, in order 
to detect the relevance of the two-point entanglement versus 
the n-way entanglement (n > 2), the "CKW conjecture" was 
assumed, i.e., the conjecture that the inequality (fT4t can be 
extended to states of an arbitrary number of qubitsi^ 

ti = C 2 A ( BC > C A b + C\c + • • ■ + C AN — r 2 . (16) 

Note that for spin systems all the concurrences can be easily 
evaluated due to the qubit nature of the subsystems. In Ref. 

starting from the above conjecture, the authors define the 
entanglement ratio Er as follows: 

E R = r 2 /n < 1 . (17) 

The more the ratio decreases, the more QS correlations are 
relevant with respect to Q2 ones. Recently the CKW conjec- 
ture has been rigorously proven in Ref. |24j. 



Generalizations of the above results were carried out in 
Refs. |26|and|27l where the authors provided a bipartite entan- 
glement measure for the case of arbitrary dimensions of the 
subsystems by defining, even in this case, the notion of tangle. 
The latter construction is, in general, difficult to apply, since 
the determination of the generalized bipartite concurrence re- 
quires application of optimization processes. In our case, this 
implies that Er can be easily applied only in the phase where 
the local Hilbert space is of a two-qubit kind, i.e., in phase I 
and III. 

In order to overcome this problem and study the transitions 
II — > I,III,IV, we will make use of a different kind of ratio that 
allows one to compare the total two-point correlations with 
the total correlations (quantum) of a single site with respect to 
the others. This can be done by substituting in the definition of 
the entanglement ratio: (i) the sum of the squares of two-site 
concurrences r 2 with the sum of quantum mutual information 

L-l 

(ii) the linear entropy t\ with the single-site entanglement 
f i = S L . The new ratio, termed correlation ratio, reads 

C R = f 2 /n.. (19) 

The fact that quantum correlations cannot be freely shared by 
many object is encoded in the so-called monogamy principle 
demonstrated in Ref. 24, which is at the base of the definition 
of r 2 . Classical correlations C2 are not required to satisfy 
this principle; hence the sum of the mutual information of a 
given site with the remaining of the lattice is, in general, not 
bounded. Such a feature, however, does not affect the change 
of Cr at QPTs. Cr compares the two-point correlations (Q2 
+ C2) of the site i with its total correlations, that, in our case, 
are purely quantum (Q2 + QS). As we shall show, it is a 
useful tool to characterize the phase transitions II — * I,III,IV 
in terms of two-point versus shared correlations. 



IV. RESULTS AND DISCUSSION 

In this section we first derive in the p — u setting the re- 
sults achieved by the method presented in Ref. 10 about the 
two-point and/or multipartite nature of the entanglement in- 
volved at each transition for the model Hamiltonian Hrc- We 
then deepen the analysis by employing the measures described 
above [Eqs. ( fTTT ). ( fTTT i and ( fT9] ll at the same transitions. 

The method described in Sec. IIHAI classifies the type of 
entanglement involved at a given QPT by direct comparison 
of the derivatives of Si and . In order to evaluate these two 
quantities, the eigenvalues of the reduced density matrices pi, 
and pij are needed. While the on-site reduced density matrix 
(O is already diagonal, the two-site density matrix (0 is block 
diagonal and, in its diagonal form, reads 

p ij =diag{D 1 ,0 + ,0-,D 2 ,V + ,V-,Q + ,Q-,D 3 }, (20) 
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where, in the thermodynamic limit (TDL), 

o± = (i-c)k(i-n,) + |r y |(|ry|±i) 



Q+ 
Q- 



2c(l 
0, 



(1 



|r*il 2 



V± = c 



n s {l-n s ) + \T ij \(\T ij \±l) 



(21) 

(22) 
(23) 
(24) 



In Table HI] we summarize the behavior of the various func- 
tionals evaluated at the transition points. In the last column, 
the transition is labeled as Q2 whenever the divergencies dis- 
played by d x Si and d x Iij are of the same type. In the other 
cases, i.e., when only d x Si displays a divergency, the transi- 



tions are labeled QS. As already discussed in Sec. IIII Al the 
two groups reflect the relevance and the role of the Q2 and QS 
correlations at the different transitions. The detailed analysis 
carried out in the following sections not only fully confirms 
the existence of these two groups, but also gives evidence of a 
further unexpected critical phenomenon occurring at each of 
the transition of Q2 type. 

A first consideration about Table [TT] is in order. As already 
mentioned, we consider the derivatives with respect to u and 
fj,, which are the quantities that parametrize the Hamiltonian 
(H). At variance with the study in the u — n setting developed 
in Ref. 10, the transition I,F — > IV is here described by d^Si 
and dfxTij. This allows us to properly include it in the Q2 
group. 















d x E R 


Ent 


1,1' -f IV (x -■ 


= M) i/VIm - 




-1/tt 2 






i/^\n - ftc\ 


Q2 


ii -> i,r (x = 


u) log(u c — u) 


f 


f 




f 




QS 


II -> III (x = 


u) l/y/u — U c 


l/\/u — U c 


l/(2n 2 ) 




l/^U — Uc 




Q2 


II -> IV (x = 


u) 1/Vmc — u 


1/V u e — U 


-1/(4tt 2 ) 




1/V u e - U 




Q2 



TABLE II: Behavior of Rjj- and the evaluated partial derivatives at the various QPTs (left column): the critical values u c and /j, c can be inferred 
from TableU "f ' stands for finite value. 



A. Two-point entanglement at Q2 transition points 

In this section, we analyze the behavior of the two-point 
correlations at each of the Q2 transitions. We first proceed 
with the computation of the negativity (fTTT) . At variance with 
the concurrence, the negativity can be used even when the lo- 
cal subsystem's Hilbert space has dimension greater than 2, 
i.e., region II. We apply the partial transposition to pij, and 
then we proceed with its diagonalization. pTj is block diago- 
nal and the only nondiagonal subblock reads 



where 7^ 



1 



(l-O ± y/1 - 2n a (l - n a ). 



Quite interestingly, the result shows that the region charac- 
terized by nonvanishing negativity, i.e., r?- G [7^,7+], de- 
pends only on n s . Such a result reveals that the presence of 
Q2 correlations in the TDL is deeply connected with the pres- 
ence of the fermionic Sutherland specie: when the latter is 
absent, all negativities go to zero, whereas when n s 7^ 0, the 
relative number of empty or doubly occupied sites does not 
influence the presence of two-point entanglement. 



2 D 2 P 2 
V Q Pi D 3 



(25) 



The negative eigenvalues of p t j coincide with the negative 
eigenvalues of this subblock. In the TDL, the only possibly 
negative eigenvalue reads 



X. 



Pij +D 2 -J [aP l3 - D 2 Y + 4a|r 2 



(26) 

with a — c 2 + (1 — c) 2 . A straightforward calculation shows 
that 



M -a_, 7 i<rj< 7i 

otherwise , 



(27) 



1. Concurrence versus Negativity 

Before proceeding in our analysis of the behavior of Q2 
correlations at the transition points, it is useful to compare 
results obtained through negativity ( 1271 with those obtained 
with concurrence (TlOb in the regions where the latter can be 
evaluated, i.e., regions I and III. As we shall see, even though 
it has been proven that both are measures of entanglement for 
qubit systems, the comparison shows that, in general, they 
have different behaviors and derivatives. As far as phase I 
is concerned, we have that the concurrence is given by 



CL = 2max 0, IF 



[(1 



|r y | 2 ] (n a -|r«p) 



(28) 
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whereas specializing Eq. (|27T i to the case of region I, we ob- 



tain 



Affj = max|o,i [(l-n) 2 + n 2 -2|r t 



(l-2n)2 + 4|r tf |2 



(29) 



Figures |2] and [3] show that, as expected, both measures are 
nonvanishing in the same intervals [see Eq. (1271) 1. Indeed, for 
r = \i — j\ > 1, they start to differ from zero in correspon- 
dence with the same values of /i. When both concurrence 
and negativity are nonzero, their behavior differs in at least 
two relevant aspects. First, apart from the case r — 1, the 
two quantities reach their maximum in correspondence with 
two different values of fi. Such a feature is compatible with 
the fact that, in general, the two measures provide a differ- 
ent ordering of the states. 22 Second, they differ in the behav- 
ior of their derivatives with respect to /j,. In particular, while 
at transition I— > IV the derivative of concurrence does dis- 
play the correct diverging behavior (d^Cij w l/yj/Xc — /x|), 
dfj,Nij does not display any divergence. As for region III, in 




FIG. 2: Region I, u = 4. The curves a and c are the concurrencies 
Ci and C2, respectively, while the curves b and d are the negativities 
Ai and A/2, respectively. 

the TDL both Afij and are always zero. The behavior of 
the two quantities significantly differs if finite-size effects are 
included. Indeed, to first order in 1/L and for all \i — j\, the 
concurrence reads 




n d ^0 
otherwise , 



(30) 



2. Divergence of the entanglement range 

The fact that A/y differs from zero at different values of \x 
depending on r = \i — j\ allows one to identify the range of 
negativity 

Rx{u,ii) = {r I Af iii+r ^ A M.i+r+i = 0} , 

i.e., the maximum distance r at which the Negativity is non- 
vanishing at fixed u, fi. We find that Rj^ is always finite ex- 
cept in the two following situations: (i) when n s — > (tran- 
sition II — > III) and (ii) when n s — > 1 (transitions I — ► IV 
and II — > IV), in which R^j- diverges. In particular, R^/ re- 
mains finite in correspondence with the two transitions II — > 
I, at which n s — * n ^ 1 [and correspondingly u — > u c (n)]. 
In these cases, only the nearest-neighbor negativity is always 
positive, i.e., Rj*j- > 1. The condition that fixes the value of 
Rtf is again Tj = 7 2 . 

The entanglement range allows one to better characterize 
the difference between Q2 and QS transitions. In the lat- 
ter, the generic site i is correlated — via two-point quantum 
correlations — only with a finite number of neighboring sites, 
whereas, at Q2 transitions, the two-point quantum correla- 
tions begin to spread along the chain and, at the critical point, 
two arbitrarily distant sites are quantum correlated. This latter 
case is shown in Fig. [3] for the transition I — > IV, at which 
Ca can be evaluated as well. We can further characterize 




0,05 .0,04 .0,03 .0,02 -0,01 0,00 



FIG. 3: Region I, u = 4. Aij (right) and Cy (left) for r = 1, .., 8. 

the spreading of the correlations by analytically exploring the 
scaling behavior of Rfj at the transitions. We have that 




(32) 



where Rn rj 0.44 is the solution of 



while the negativity i 



;„10 



yw-III = n d {l-n d ) 1 
n 2 + (1 - n d y L 



sin 7ri?o 1 



(33) 



(31) The exponents v c characterizing the divergence of Rj^ at u c 
and fi c for the various transitions are easily worked out from 
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Eq. d32b by recalling that n s is a function of u and /x (see 
Table Q). Quite interestingly, at all the three Q2 transitions, 
we have 



(34) 



A similar type of behavior was already studied for a spin 
model 7 where the notion of entanglement transition was in- 
troduced. In that case, the divergence of the entanglement 
range is observed for the concurrence at some specific point of 
the phase diagram for which the ground state becomes factor- 
ized, and apparently no QPT takes place. We recognize such 
a feature for the model discussed here only in part. In fact, 
here all the entanglement transitions occur in correspondence 
with QPTs; moreover, while phase IV is indeed characterized 
by a factorized structure (the ground state being singly occu- 
pied at each site), phase III is not, since the ground state in 
this phase is a superposition of empty and doubly occupied 
states distributed over the whole chain. This observation sug- 
gests the conjecture that a factorized structure with respect to 
Hilbert space appearing in the ground state is a sufficient but 
not necessary condition for the occurrence of an entanglement 
transition. The conjecture could also be generalized in terms 
of bipartite entanglement: an entanglement transition occurs 
if and only if the new phase is a two-point entanglement free 
one. In this sense the factorized state of phase IV and the 
genuine multipartite ground state of phase III are equivalent. 
Moreover, at least in our model, it is equivalent the way in 
which the system destroys all correlations (IV) or build gen- 
uine multipartite ones (III). 



B. Two-point versus multipartite entanglement at QPTs 

In order to explore the role of multipartite entanglement at 
the various transitions, an ideal tool would be the entangle- 
ment ratio Er d 17T >. which, as explained, provides a direct 
measure of the relative role of Q2 and QS: a decreasing (in- 
creasing) entanglement ratio in proximity of a QPT means that 
QS (Q2) correlations are more relevant to the transition. Ac- 
cording to Eq. ( TToT l. T2 is properly measured through con- 
currence. This implies that only the transitions in which the 
system is of qubit nature (namely, the I— > IV transition) can 
be explored through Er, whereas in region II we used the 
correlation ratio. 



requires some attention. In fact, Cy depends on the distance r 
between the two sites; one has to first evaluate the sum 



L-l 



r 2 (iV s ,L) = ^C 2 (7V s ,L) 



(36) 



and then the TDL limjv 3 ,L-«x> T 2 {N s , L). The numerical eval- 
uation of En for a sufficiently large r is reported in Fig. [4] 
The latter clearly shows that, as expected, in the vicinity of 




FIG. 4: Region I, u = 4. Entanglement ratio Er and normalized 
correlation ratio Cr/L with L — 1000. 

transition I — > IV (p — > 0) the Q2 correlations rapidly in- 
crease with respect to QS ones. We notice that interestingly 
the derivative of Er diverges again with |/i| -1 / 2 at the transi- 
tion. 

We pass to explore r 2 in region III, which is again of qubit 
nature. Since such a region is characterized by the presence 
of just j] pairs, which have an intrinsic multipartite nature, we 
expect Er to vanish there. As far as t\ is concerned, we have 
that 



Til 



4 det [p 



IIIi 



4n d (l - n d ) 



(37) 

whereas for r 2 one can see that, since 

4 IIT) ~ l/L [see Eq. 
(f30b l is independent of r = \i — j\ and vanishing in the TDL, 
t 2 = J2j C?j — 0- Hence, the entanglement ratio correctly 
indicates that the only relevant entanglement is multipartite 
(i.e., rt-way entanglement with n > 3). 



1. Entanglement ratio 

We start exploring the behavior of Er in phase I, where Cy 
is defined, and, in particular, at the transition I — » IV. As far 
as n is concerned, from Eq. ( fTBT l we have 

t[ = 4 det \p{ ] = 4n(l - n) . (35) 

As for T2, the two-site concurrence was given in Eq. d28l ); 
still, the evaluation of the sum of the {L — 1) Ch in Eq. (fT~6b 



2. QS correlations in region II 

As a general fact, in region II both Q2 and QS correlations 
are present. In particular, the QS transitions — which accord- 
ing to Table HIl are II — > I,F at fixed u — should be character- 
ized by some change in multipartite entanglement QS. Such 
a hypothesis has a first strong confirmation in the fact that at 
these transitions 77 pairs [and ODLRO, see Eq. ([5])] disappear. 
Indeed, it has been shown in Ref . l3l| that r\ pairs do carry mul- 
tipartite entanglement, thus disappearing at these transitions. 
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Furthermore, the behavior of Q2 correlations is here radically 
different from the one they display at Q2 transitions. Actually, 
as already seen in Sec. IIV A 21 the entanglement range is finite 
at any filling n < 1 for u < u c (n); moreover, both and 
Mij, and their derivatives remain finite in the same regime. 
At n = 1 and u = 4 (i.e., at transition II—* IV), besides QS 
also Q2 correlations enter to play a role. In fact, while r\ pairs 
disappear, Rj^ becomes infinite and an entanglement transi- 
tion takes place. Since the analysis of the previous section has 
shown that R^/ has the same divergence of Si, we infer that 
the role of Q2 correlations is dominant at this transition. In 
order to confirm such a scheme, we now use the correlation 
ratio previously introduced. 

3. Correlation ratio 

We aim at obtaining an indicator of the relative weight of 
Q2 correlations with respect to QS ones in region EL Since Zy 
keeps track of the change of Q2 correlations between i and j 
at transition points, we expect the correlation ratio Cr ( fT9l ) to 
capture such a desired feature. 

We first consider the behavior of Cr in region I, where it 
can be compared with the standard entanglement ratio Er; 
this is shown for sufficiently large L by the dashed line in 
Fig. [4] In correspondence with the transition II — > IV, Cr 
correctly reproduces the qualitative behavior of Er, i.e., the 
relative weight of Q2 correlations rapidly increases. In Figs. 




FIG. 5: Region II, n = 1/2 u G [—4,0]. Normalized correlation 
ratio Cr/ L with L — 1000. Inset: zoom of the transition II — > I. 

|5]and|6]we report Cr(u) in region II at two different values of 
n (n — 1/2 and n — 1). At both values Cr rapidly increases 
in proximity of u — — 4 (transition II — > III). The behavior of 
Cr is quite different in the two cases in correspondence with 
the upper critical point. Indeed, for n = 1/2 (transition II — > 
I, Fig. it goes to zero with a clear linear dependence on 
u c — u, reminiscent of the behavior of the pairing gap Ap, 
whereas for n = 1 (transition II — > IV), after decreasing in 
almost the whole region, Cr rapidly increases for u — > 0. 
These features are in accordance with the considerations 



0,20- 



0,15- 



0,10- 



0,05- 



0,00 




u 



FIG. 6: Region II, n = 1, u 6 [—4, 4]. Normalized correlation ratio 
Cr/ L with L = 1000. Inset: zoom of the transition II — > IV. 

exposed in the previous section. Whenever Cr increases at 
the transition, this implies that Q2 correlations are increasing 
with respect to QS, and hence the transition should be of Q2 
type. On the other hand, it is only at the transition of QS 
type (transition II — > I) that Cr vanishes, meaning that QS 
correlations overcome Q2 ones. 



C. Entanglement away from QPTs 

Apart from transition points, we can spot the areas in the 
different regions where the QS correlations prevail with re- 
spect to Q2 correlations from the direct study of Si, Mij. 
In Region III, QS correlations prevail everywhere since Si 
is different from zero and all A/y are vanishing. Let us re- 
call that one has genuine multipartite entanglement whenever 
both Si 7^ and two-point entanglement is zero. Here it is 
due to the presence of r\ pairs, which is also captured by two- 
point classical correlations. In fact, it turns out that l\j ^ = 

l£ n) = 2n d (l - rid) Vj, with Zoo = lim| i _ 7 -|_ ) . 0o I\i-j\. 
All pairs of sites are equally correlated as two infinitely dis- 
tant sites. Interestingly this property is directly related to the 
presence of ODLRO in that the total amount of correlations is 
simply proportional to it. 

In region I, QS correlations prevail away from transi- 
tion points since the entanglement ratio has a minimum (see 
Fig. |4|i. Contextually, only the nearest-neighbor negativity is 

nonzero, and Si is maximum; moreover, Zio = 0. The same 
qualitative behavior holds inside region II as well, except that 
loo 7 ' 7^ in the whole region except at the transition II — * I, 
as can be seen from the dashed line in Fig. [8] This is related 
with the fact that rj pairs are present in region II as well. Quite 
interestingly, the contribution of singly occupied and doubly 
occupied sites to two-point correlations seems to simply add 
in quantum mutual information. Indeed, one could check that 
T (H) _ T (j) T (H) 

To resume, we observe that an infinite range of two-point 
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correlations in proximity of a transition is a signal of a Q2- 
driven QPT, in which case R^j- also diverges. Far from transi- 
tion, the same infinite range is implied whenever ^ and 
it is thus a signal of the existing ODLRO, in our case related 
to 77 pairs. 




FIG. 7: Region I, u — 4. Mutual information , r = 1, .., 5. 



1 1 1-3 I 
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FIG. 8: Region II, n = lit 6 [—4,4]. Mutual information 2, 
r = l,..,5. 



V. CONCLUSIONS 

In this paper, we analyzed the rich phase diagram of the 
one-dimensional bond-charge extended Hubbard model at 
T = by means of various measures of bipartite and/or mul- 
tipartite correlations. All the computed measures are capable 
of reproducing the known phase diagram in terms of singu- 
larities; moreover, at each transition the critical exponent of 
the correlation length is shown to coincide with the scaling 
exponent of the divergent quantities, when evaluated. 

The knowledge of one- and two-site Von Neumann en- 
tropies allows one to distinguish the quantum phase transi- 



tions (QPTs), in which the role of two-point quantum correla- 
tions (Q2) is relevant to those in which multipartite quantum 
correlations (QS) are determinant. 

The systematic analysis of the appropriate measures of en- 
tanglement at the different transitions and phases made it pos- 
sible to better characterize both Q2 and QS transitions. Con- 
textually, a different estimator which can be computed for qu- 
dit systems has been introduced: the correlation ratio. The 
latter is an indicator of the relative weight of the two-point 
correlations with respect to the total (quantum) ones at the 
transitions. The analysis shows that Q2 transitions are char- 
acterized by a divergence in the range of negativity: when 
approaching the transition, two sites at arbitrary distance be- 
come quantum correlated. Such a feature is reflected in the 
behavior of the entanglement and/or correlation ratio. The 
latter increases at Q2 critical points with diverging derivative, 
clearly indicating an increasing relative weight of two-point 
quantum correlations with respect to multipartite ones. QS 
transitions instead are characterized by a finite range of nega- 
tivity and by a vanishing correlation ratio, indicating that mul- 
tipartite quantum correlations dominate there. For our model, 
the correlated physical phenomenon is the disappearance of 77 
pairs. 

Finally, we described the nature of the correlations within 
each region as well. For our model, the existence of two-point 
quantum correlations depends only on the presence of singly 
occupied sites. At the same time, the presence of doubly oc- 
cupied sites witnesses the appearance of 77 pairs and ODLRO, 
and multipartite entanglement carried by them. At the level 
of two-point correlations ODLRO coincides with the finite 
value of quantum mutual information between infinitely dis- 
tant sites. 

In conclusion, the above analysis has widely clarified how 
to characterize the nature of quantum correlations involved 
at a QPT for an integrable correlated electron model. The 
scheme, in particular, allows one to gain from quantum mutual 
information insight on the behavior of both Q2 and QS corre- 
lations at transition points for T = 0. We expect the scheme 
to be straightforwardly applicable also in nonintegrable cases, 
both in one and in greater dimension. A first step in this 
direction has been achieved in one dimension by means of 
the numerical analysis in the nonintegrable case ([]]). — It re- 
mains to be investigated how to modify the proposed scheme 
at T 7^ 0, where also temperature-driven correlations play a 
major role; in particular, it is expected that they would com- 
pete with quantum ones in determining the behavior of quan- 
tum mutual information. 
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APPENDIX A: REDUCED DENSITY MATRIX EXPANSION 



a. Region I 



Let us recall the definition of the reduced density matrix 
PU- 

Pij=TrL/i j \*Gs)(*Gs\; (Al) 
given that the pure state |^gs) can be written as 

|* GS ) = ^2Cab\a)\b), C ab eC, (A2) 

ab 

where {\a}} and {\b)} are the basis for the subsystems {ij} 
and {L/ij}, respectively, the operator expansion for pij is 
easily derived: 



Pij = T^t/H E C ab C* a , b ,\a)\b){b'\{a'\ 

aa'bb' 

= E C ab C* a , b/ (b"\a)\b)(b'\(a'\b") 

aa'bb'b" 

= J2 C ab C:, b ,(-)( ab+a ' b ')\a)(a'\5 bb// 6 b , b „ 

aa'bb'b" 

= J2 c ^c:> b (-)^ +a >\a)(a'\ 

aa'b 

= E c -»( a 'i » < A3 > 



where 



GS 



(A4) 



where (— ) a takes into account the parity of the state \a)M- It 
turns out that for our model a + a' = for all (a, a'), due to 
the conservation of N s and Nd- 



APPENDIX B: DIMER REDUCED DENSITY MATRIX 
EVALUATION 

In this section, we schematically give the procedures to 
compute the dimer reduced density matrix. The mean value 
of the following operators give the diagonal elements of p^j : 



y02 v20 y02 y20 y02y20yll Ylly02y20 

A i A i A j A j ' A i A i A j ! A i A j A j ! 

Yllyll v\1v22 y22 yll 

A i A j ) A i A j 7 A i A i ' 

V02 v20 v22 v22 v02 v20 y22 v22 

A i A i A . ) A ( A . A , 7 A i A j • 



(Bl) 



Due to the conservation of N s and Nd, the only off-diagonal 
elements that can be nonvanishing in some of the regions are 
the following ones (together with their Hermitian conjugates): 



\^10^01 y21 v 12 y02 y20 

A i A j ) ~ A i A j ) A i A j • 



(B2) 



The action of the latter is simply to permute the states on the 
two sites. 



Since the ground state in region I is given by a superposition 
of states in which each site is empty or singly occupied and 
since N s = N tot = N and XfXf = [1 - X} 1 ] the only 
nonzero entries are given by 

{xfxrxfx?^ ((i-xm-xj^, 

(X™X?°X}\ ee ((1 - X^)X}% 
(X^XfXf^ ee {X}\l~Xf)) h (B3) 



and 



(X^X}% (X}°X°\ ; (B4) 



where (O)t ee ($i(N, L)\0\^i(N, L)). We thus compute 
(Xl 1 ^ using the expression of both | * j (iV, L)) and X} 1 in 
terms of momentum operators 

< X ">' = I E ex P« fc " k , )j}(Xl°X^) I = ^ , (B5) 



k.k> 



are 



since the only nonvanishing terms of the sum in Eq. 
those for which k — k! = 0. 

The calculation of (X}° X® 1 )^ is analogous to the previous 
one; here we have the appearance of a phase factor: 

(Xl X°\ = ^ E ex PH2?rm(i - j)/L) = . (B6) 

m 

Such a phase factor has a different expression for N even or 
odd: 



|rf_,| = |r?_,.| = sin( ^ 



ir\i - j\) 



Lsin(^n\i- j\) 



(B7) 



In the TDL, = Tfj = and it can be computed by ap- 



proximating the sum in Eq. ( IB6b with the following integral: 

/>7717r 



1 



7T 



sin (n\i — j\Tr) 
cosk{i- j)dk = y ^ L - L . (B8) 







TT\t-J\ 



The previous calculations allow us to eventually compute 

(X?X]\= ^ -r^Ti-i. (B9) 

The expressions of D\, D2, 0\, and O2 in region I (Nd = 0) 
follow from the collection of the previous results. 



b. Region III 

The ground state in region III is given by a superposition 
of states in which each site is empty or doubly occupied. The 
only nonzero entries are given by 



/ V Q2 v 20 v 02 v 20\ 
( A i A i X j X j 



(X° 2 Xf°X 22 ^ 



/IIL 



{xfxfxf\ 



III' 



(BIO) 
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and 



/ v 22^22\ / v 02 v 20\ 

\ X i X j )llh \ X i X j /III > 



(Bll) 



where ( Q) m = (^ m {N, L)\0\^ m (N, L)). The evalua tion 
of Eq. dBlOt follows from the evaluation of Eq. (IB1 H . A 
general strategy is used. In the case of (X 22 X 22 ) jjj one has to 
count the number of states in the superposition ^^(./V, L)) 
whose sites i and j are doubly occupied: 



(X?X?) m = N m (N d ) 



L-2\ = N d {N d - 1) 
N d - 2) L(L - 1) 



(B12) 



In order to compute (X® 2 X 2Q )y\\ one has to count all the 
states whose site i is doubly occupied and whose site j is 
empty. This leads to the following result: 



( X T X ? >III = %l(^) 



L - 2 
N d -1 



N d (L - N d , 



L(L-l) 



(B13) 



The expressions of D\, -D3, and Q in region III {N s = 0) 
follow from the previous results. 



c. Region II 

We start by noting that since [Xf 1 ,^'] = 0, the operator 
X} 1 does not affect the doubly occupied part of the ground 
state. Accordingly, 

(N A \X?\N d ) u = ^, 



(A^xflA^n - -ir '- 



1 ' 



where the notation {N d \0\N d )n = 
(* n (JV a ,JV l i ) L)|0|* n (iV 5 ,iV < ,,L)) will be useful in 
the following calculations. 

We now compute ( N d \ X 22 1 N d ) . In a first step we make use 
of the following relations: 

(i) [Xf t ,tf) lf ']=N i tf)»*- 1 X*>, 

(ii) [Xf, (r,t)^-i] = 0, 

(iii) ( V ) N * = {rj) N *-\ 
The latter imply that 

[N^N^i^Kv^XfivY^i) 

= [N ll (N d )] 2 N d (* l \( v ) N «tf) N «- 1 Xf\y l ) 



^02 y2Q 

may write the following recursive equation: 



If we now define £V d _i = (N d - l\X° 2 X 20 \N d - 1), we 



DNj — 1 ; 



N d 



— Djvw-i, (B14) 
L L-N s -N d + 1 d ' 

whose solution is V^^ = . Thus, by collecting the 

[AW)] 2 



'N d 

above results, we have 



(N d \X 22 \N d ) = 



N, 



[Nn(N d - 1)] 



iNePm-x = -r (B15) 



We now compute (N^X^ 1 X 22 \N d ) by resorting to the solu- 
tion of a recursive equation. Since [X 22 , X} 1 ] — 0, we can 
first apply the procedure used for (N d \X 22 \N d ) to obtain 



(NdlXfX^lNd) 



Wn(N d ) 



[N n (N d -l)}' 



:(N d ) 2 £ Nd - U (B16) 



^02 Y 20 y11 

following recursive expression 



where £ Ni _ x = (N d - l\Xf Xf Xj^Na - 1) satisfies the 



£ Nd -i = (Ni-mi-Xj^XflNa-l) 



[AW -I)] 2 
[N n (N d -2)} 2 

N N 2 

— - — i + ir- ■ 

N d -1 



(N d - l) 2 £ Nd -2 



L-N s -N d + 2 
The solution of the latter is 

N d 



-£ 



N d -2 



'N d 



1 



L-N s 
and, finally, we get 



(NdlXfX^Nd) = 



L 



L 



(B17) 



01 8) 



N d 



L - N* 



N S (L-N S ) 
L 2 



- »— j 1 



B.19) 



We now compute (N d \X 22 X 22 \N d ). The same arguments 
used for Eq. ( IB 17b lead to 



(N d \X 22 X 22 \N d 



[AW)] 2 
[%(^-2)] : 



where the function defined as Tm, 

{N d \X° 2 X 20 X° 2 X 20 \N d ) satisfies the followil 
sive equation: 



T m = (a - a') ( 1 - 



P J P-m+l 



T m -x , (B21) 



with 



a = 1-^, p = L-N s 



f -7^ + ir^l 2 , - 



(B22) 



The latter recursive equation is solved by defining the auxil- 
iary function £ m = <3 _" +1 ^ r m -i ( a i a >) that obeys 

1 , (B23) 



p-m+r 



whose solution is G m = — ■ Collecting the above 

results, we have that 

(N d \X 22 X 22 \N d ) 



(B24) 



N d {N d - 1) 


;(i-f) 2 -iiVii 2 ; 


(L-I\ 


/.)(i-JV.-l) 
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The computation of (N d \X° 2 X 20 \N d ) is now straightfor- 
ward, since it can be expressed in terms of the above defined 



(N d \X° 2 X 20 \N d ) = 



N d (L — N s — 


N d ) 


> 


-f) 2 -|r^-| 2 ; 


(L- 


N S )(L- 


N.-l) 



(B25) 



The same argument holds for (N d \XfX 20 Xf\N d ): 



{N d \XfX™X?\N d ) = 



W^^W d Nd -' 



= (N d \X° 2 X 20 \N d 



(B26) 



The previous steps allow us to easily evaluate 
(N^XfX^X] 1 \N d ) in terms of Eq. (iBTgl : 



(N d \X° 2 X 20 X} 1 \N d )=£ Nd . 



(B27) 



We then compute (NdlX^X^lNd). Using the following re- 
lations: 



(i) [X^Xf , fat)"-] = -N d X}°Xf (rf) N *-\ 



(ii) {rf) Nd X 21 = N d (r]) Nd ^ 1 X° 1 ; we obtain the recursive 



equation 



£(JV d ) = (N d \Xl°X^\N d ) 
N d C(N d - 1) 



N d - L - N s + 1 



+ IV, , (B28) 



whose solution is 



{N d \Xl X° 1 \N d )=T 



i-3 



' 111 s ' nli N d ) = r L ~ N <i 

-12 y-211 



i — iV. 



(B29) 



The task of evaluating (N d \X} 2 X 21 \N d ) is simplified by ob- 
serving that [XfX^rf] = [X? 1 Xj 2 , 7jt] ; we obtain 



r^(iV d - 1) 



N,, 



-r, 



L-N s ' 3 -'" 
where C(N d — 1) is given by Eq. dB281 >. 



(B30) 
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